Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 10 May 1995 



Extensions around z = 2 QSOs 

Itziar Aretxaga, B.J. Boyle, Roberto J. Terlevich 

Royal Greenwich Observatory, Madingley Road, Cambridge CBS OEZ, U.K. 



Submitted 30 March 1995 



On 
On 



ABSTRACT 

Based on an R-band imaging survey of 6 high redshift (z ~ 2) and high luminosity 
(M B ^ —28 mag) QSOs, we report the detection of extensions in two radio-quiet and 
one radio-loud QSO. The extensions are most likely due to the host galaxies of these 
QSOs, with luminosities of at least 3 — 7 per cent of the QSO luminosity. The most 
likely values for the luminosity of the host galaxies lie in the range 6 — 18 per cent of 
the QSO luminosity. 
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©1. INTRODUCTION 
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The study of high redshift (z ~ 2) QSOs offers a unique 
opportunity to investigate conditions in the early universe. 
In the currently favoured cold dark matter cosmogony, 
this epoch corresponds to the period when normal galax- 
ly-^ ies formed through hierarchical coalescence (Carlberg 1990), 
Q thereby giving rise to enormous concentrations of gas in the 
\f*} center of the galaxies, which could feed a central black hole 
(Rees 1988) or provoke a massive starburst episode (Ter- 
levich & Boyle 1993). As such, this picture is consistent with 
the observation that QSOs exhibit a peak in their activity 
I at these redshifts (Schmidt, Schneider & Gunn 1991). In- 
O deed, there has been much observational effort to obtain 
+3 deep multicolor images of high redshift radio galaxies and 
^radio-loud QSOs with spectacular results: enormous struc- 
^ tures of up to 100 kpc aligned with the radio-jet (Chambers 
et al. 1987), luminosities and sizes of the parent galaxies 
larger than those of normal elliptical galaxies (Lehnert et 
al. 1992). However, radio-loud quasars are only a small frac- 
tion (< 5%) of all QSOs, and many of the conclusions drawn 
from the studies of radio-loud objects might be unrepresen- 
tative of the conditions in the early universe. It is, therefore, 
important to look into the properties of radio-quiet QSOs, 
as a better indicator of the general properties of galaxies 
at high redshifts. We present here deep images of a sample 
of 6 high redshift, predominantly radio-quiet QSOs. Section 
2 gives a summary of the observations, section 3 explains 
the technique used to detect extensions around the QSOs, 
and section 4 discusses the likelihood that the extensions are 
host galaxies. 



OBSERVATIONS 



band with the f/11 auxiliary port of the 4.2m William Her- 
schel Telescope (WHT) of the Observatorio de Roque de los 
Muchachos in La Palma. Luminous (Mb < — 28 mag)* high 
redshift (1.8 < z < 2.2) QSOs which lay close (20 £ 6 <> 
60 arcsec) to stars of similar brightness, were selected from 
the Veron-Cetty and Veron Catalogue (1994) catalogue. The 
existence of a nearby star allowed us to determine an accu- 
rate point spread function (PSF) for each CCD frame. The 
observations were performed under excellent seeing (< 0".9) 
in dusty and smoky non-photometric conditions. 

We carried out the observations with a 1024 2 Tektronix 
CCD chip (TEK1), giving a spatial resolution of 0.105 arc- 
sec/pixel over an unvignetted 1.8 arcmin diameter field. The 
configuration was chosen to provide an optimal sampling of 
the PSF, even in the best seeing conditions (« 0".5). The 
chip was readout in QUICK mode, which yields a read noise 
of 4.3 e" and a gain of 1.8 e"/ADU. 

Each QSO field was observed in a series of 600-900s ex- 
posures to avoid saturation of the QSOs and PSF stars. Indi- 
vidual exposures were offset from one another by ~ 10 arcsec 
to prevent the repetition of any coherent signal (sources, 
cosmic rays, bad pixels) at a particular location on the chip. 
Table 1 summarizes the observations. 



3. DATA REDUCTION AND ANALYSIS 

The data were reduced using the IRAF software package. 
For each separate QSO field, the images were first bias sub- 
tracted and then flat-field corrected using the sky flat-field 
derived from the median of all individual frames, once bright 
stars, galaxies and QSOs (and their surrounding regions) 
had been masked from the original data. 



Four radio-quiet and two radio-loud QSOs were observed 
during 1994 August 8/9 and 9/10 in the Harris R pass- 



* Ho = 50 Km s 1 Mpc 1 and qo = 0.5 are used throughout 
this work. 
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Table 1. Summary of observations 



Name 


B* 


M B * 


z* 


Radio-loud?* 


Exposure time 


No. of frames 


FWHMt 












(seconds) 




(arcsec) 


UM 281 


17.4 


-28.9 


1.878 


N 


6600 


11 


0.7 


1630.5+3749 


18.3 


-28.3 


2.037 


N 


10500 


15 


0.6 


PKS 2044-16 


17.4 


-29.1 


1.943 


Y 


5400 


9 


0.9 


PKS 2134+008 


16.8 


-29.6 


1.936 


Y 


6300 


7 


0.6 


Q 2230+02324 


18.0 


-28.7 


2.150 


N 


5400 


9 


0.6 


Q 2244-0105 


18.0 


-28.6 


2.040 


N 


5400 


9 


0.7 



* B, Mb, z and radio properties from Veron-Cetty and Veron (1994). 
t Stellar FWHM measured on co-added frame. 



We then applied the PSF subtraction techniques pro- 
vided by the IRAF package DAOPHOT to establish whether 
QSOs in this study exhibited any significant non-stellar ex- 
tensions to their profiles. This procedure has already been 
successfully used to isolate intervening galaxies in the line of 
sight of QSOs (Steidel & Hamilton 1992). For each individ- 
ual frame, we defined a 2-dimensional PSF using the bright- 
est of the closest stellar companion to the QSO. The PSF 
was defined over 30-40 pixels in radius. We then scaled the 
PSF to match the luminosity of the QSO and other nearby 
stars over the same region, and subtracted the scaled PSF 
from them. The remaining residuals in the non-PSF stars 
provided an accurate check of the validity of the subtrac- 
tion process. We accepted the PSF subtraction if the resid- 
uals in the non-PSF stars accounted for less than la of the 
Poisson noise expected from the subtraction technique, and 
showed no systematically repeatable pattern from frame-to- 
frame. The la level was calculated by combining the Pois- 
son noise from the different sources (sky subtraction, PSF 
fitting, QSO profile) in the PSF subtraction process. For one 
field, Q2230+0232, there was just one bright star near the 
QSO which was used to model the PSF. Although the resid- 
uals for the non-PSF stars are consistent with zero, their 
distance from the field centre (> 40 arcsec) results in a dis- 
tinctive pattern in their residuals which repeats from frame 
to frame. For PKS 2044—168 there were two bright stars 
near the QSO, but neither star yielded a PSF which pro- 
duced a non-zero residual when subtracted from the other 
and, therefore, we cannot assess the significance of the resid- 
uals left in the QSO position. 

For the remaining fields the subtraction of non-PSF 
stars was successful. The subtracted frames were registered 
using the relative frame-to-frame X,Y shift derived from the 
positions of three bright stars in each original frame, and 
co-added to produce a final subtracted image. The same 
transformations were used to register and co-add the origi- 
nal frames to produce a final un-subtracted image. 

The left panels in figure 1 show the fields of the QSOs 
UM 281, 1630.5+3749, PKS 2134+008 and Q2244-0105. 
The right panels in figure 1 show the same fields after sub- 
tracting the corresponding PSFs from both QSOs and stars. 
Figure 2 shows in bigger detail 6" x 6" sections of the residu- 
als left in the positions of the QSOs and non-PSF stars. The 
residuals left in the position of the stars are consistent with 
zero at a la level (see above). Table 2 summarizes the sig- 
nificance of the residuals in the number of ADU. We have 
detected residuals in excess of 3a in the following QSOs: 
1630.5+3749 (4a), PKS 2134+008 (3a) and Q 2244-0105 



(3.7a). All the residuals show a 'doughnut' shape with a well 
of negative counts in the centre. This indicates that there 
is a flatter component below the PSF in the centres of the 
QSOs, from which the nuclear (PSF) contribution has been 
over-subtracted. To estimate the true luminosity of these 
systems, we have subtracted smaller amounts of the PSF in 
order either a) to produce zero counts in the center of the 
residuals or b) to achieve a flat-top profile with no depres- 
sion in the center. We regard these quantities as a lower limit 
and a best estimate, respectively, of the total luminosity of 
these extended components. In the lower limit case, the lu- 
minosities of the extensions comprise 3-7 per cent of the 
luminosity of the QSOs. In all cases, the FWHM of the flat- 
top residual profiles are significantly larger than the FWHM 
of the stars in each field (see tables 1 and 3). The estimates 
of the extended component luminosities derived from flat- 
top profiles increase to 6-18 per cent of the luminosity of 
the QSOs. To obtain a crude zero-point (±0.3 mag) for the 
images, we used the E magnitudes of the bright stars in the 
field provided in the APM Northern Sky Catalogue (Irwin, 
Maddox & McMahon 1994), which correspond to R band 
(Evans 1989). Since we do not know the spectral energy dis- 
tribution of the extensions, we have chosen not to correct 
for spectral slope or bandpass in these magnitudes, so that 
they correspond to magnitudes in a rest-frame band centred 
at ~A2300 A. We list these magnitudes in table 3. 

To check the reliability of the luminosities obtained 
from the PSF subtraction technique, we tried to recover 
the luminosities of known faint compact (FWHM~ 1 — 
2 arcsecond) galaxies in each CCD frame after we had 
'added in' scaled versions of the stellar PSF, mimicking 
the contribution from the unresolved bright QSO nucleus. 
We found that subtractions which produced flat-topped 
residuals could recover the input luminosity of the galax- 
ies within an error of less than 4 per cent, if the galaxies 
contributed more than 3 per cent to the total luminosity of 
the galaxy+PSF system. For galaxy contributions of & 1 per 
cent, the input luminosity of the galaxies were poorly recov- 
ered, with errors of up to 30 per cent. 



4. DISCUSSION 

Out of a sample of one radio-loud and three radio-quiet 
QSOs with suitable PSF stars, we have been able to detect 
extensions in three cases. The best estimates for luminosi- 
ties of these systems lie between 6-18 per cent of the total 
QSO luminosity. Host galaxies of radio-loud QSOs at z ~ 2 
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Figure 1. 50" x 50" fields of the QSOs UM 281, 1630.5+3749, PKS 2134+008 and Q 2244-0105. Left panels show the fields before 
subtracting the PSF stars. Right panels show the residuals left in the position of the QSOs and stars, after a luminosity scaled PSF 
has been subtracted. QSO residuals show a 'doughnut' profile, with negative counts in the center. The minimum contour level has 
been set to be ~ 2a above the sky level. North is left, and East is down. 



Table 2. Significance of the 
number of ADU 



'doughnut' residuals: 



Field 

UM 281 
1630.5+3749 
PKS 2134+008 
Q 2244-0105 



QSO 

2123 + 887 
3932 ± 1003 
7940 ± 2539 
4997 ± 1370 



PSF star 

830 ± 890 
491 ± 999 
-1491 ± 2207 
164 ± 1264 



non-PSF stars 

740 ± 754 
522 ± 733 
1706 ± 1791 
-324 ± 1071 



which comprise up to 20 per cent of the QSO R band lumi- 
nosity have already been detected by Lehnert et al (1992) 
but, to date, no other detection has been made of hosts of 
high-redshift radio-quiet QSOs. 

For radio-loud QSOs, the red colours of the extensions 



(Lehnert et al 1992) favour a host galaxy origin of the light, 
rather than scattering by dust or electrons in the halo of the 
QSOs. Although we lack colour information, the similaril- 
ity between the nuclear to extended component luminosity 
ratio in radio-loud QSOs with that observed here, leads us 
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Figurel: (continued) 



Table 3. Magnitudes of the QSOs and their extensions 



Name 


M R * (QSO) 


Mr 1 


R 1 


ADU 1 






Mr 2 


R 2 


ADU 2 






FWHM 2 


























(arcsec) 


1630.5+3749 


-28.7 


-24.8 


21.7 


15170 


± 


1001 


-25.9 


20.9 


32279 


± 


998 


1.05 


PKS 2134+008 


-30.1 


-25.9 


20.4 


45509 


± 


2536 


-26.6 


19.8 


80597 


± 


2534 


0.80 


Q 2244-0105 


-29.0 


-25.6 


20.9 


22595 


± 


1368 


-26.7 


19.9 


62975 


± 


1365 


0.84 



* Total QSO Mr absolute magnitude, including extension. 

1 Properties for extension with zero counts in the center after the PSF subtraction. 

2 Properties for extension with flat-topped profile after the PSF subtraction. 



to believe that the extensions we are detecting are also due 
to light from the host galaxy. Indeed, scattering models re- 
quire the presence of a powerful transverse radio-jet (Fabian 
1989) which is unlikely to be present in either the radio-quiet 
QSOs, or the core-dominated (Murphy, Browne & Perley 



1993) radio-loud QSO around which we have detected ex- 
tensions. 

An alternative origin for the extensions could be ex- 
tended nebular continuum, seen to be a major contribution 
in three powerful radio galaxies (Dickson et al 1995). How- 
ever, for nebular luminosities of Mr ~ —26.5 mag, as im- 
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Figure 2. 6" x 6" fields of the residuals of QSOs and non-PSF stars shown in figure 1. North is to the left and East is down. 
The centroids of the corresponding non-subtracted images lie in the center of the plots. Note that the non-PSF star in the field of 
Q 2244-0105 has a faint galaxy 1.7 arcsec to the southwest. 
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plied by these observations, the predicted narrow H/3 lumi- 
nosity would be about 3 x 10 44 erg s _1 , nearly half of the 
luminosity derived to be emitted by the broad component 
from the PSF light. Under these circumstances, the QSOs 
would exhibit prominent narrow lines with central peak in- 
tensities of more than 3 times those of the broad lines. The 
spectra of these QSOs do not show such prominent narrow 
lines (Foltz et al. 1989, Osmer et al. 1994). 

Galaxies as luminous as the extensions detected here 
have already been found in the imaging survey of radio- 
loud QSOs carried out by Lehnert et al. (1992). Four of 
the objects of their sample, with similar redshift to those in 
our sample, show 'fuzz' around the PSF of the nucleus. In 
the observed R frame the absolute magnitude of this 'fuzz' 
ranges from —24.8 to —26.3 mag, as derived from the B 
and K colours they report. Their galaxies are thus ~ 12 
times more luminous than present-day giant ellipticals and 
~ 2.5 — 3 mag brighter than the hosts of low redshift QSOs. 
Our observations show that, if the extensions we have de- 
tected are indeed galaxies, extraordinary massive and lu- 
minous galaxies are not only a characteristic of radio-loud 
objects, but of QSOs as an entire class. 

Indeed, the radio-loud QSO studied in this sample does 
not exhibit a significantly larger or a more luminous ex- 
tension than those of radio-quiet QSOs. One of the radio- 
quiet QSOs exhibits no significant evidence for any exten- 
sion. This is in contrast to the work of Lowenthal et al. 
(1995) who have established upper limits to the observed K 
luminosity of the host galaxies of a sample of high redshift 
(z « 2.5) radio-quiet QSOs to be ~ 2 mag less luminous 
than those of radio-loud QSOs (Lehnert et al. 1992). A sim- 
ilar result is found in optical studies of low-redshift QSOs 
(see Smith et al. 1986), although this difference is not ap- 
parent in the near-infrared (Dunlop et al. 1993). 

Clearly, our sample is too small to draw a firm con- 
clusion and further observations will be required to estab- 
lish the differences between radio-loud and radio-quiet host 
galaxies at high redshift. In particular, colour information 
will be vital, not only to establish conclusively that the ex- 
tensions are due to the host galaxies, but also to provide 
some clues as to the evolutionary status of such systems. 
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